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The bulk and boundary magnetizations are calculated for the critical 
Ising model on a randomly triangulated disk in the presence of a boundary 
magnetic field h. In the continuum limit this model corresponds to a c = 1/2 
conformal field theory coupled to 2D quantum gravity, with a boundary term 
breaking conformal invariance. It is found that as h increases, the average 
magnetization of a bulk spin decreases, an effect that is explained in terms 
of fluctuations of the geometry. By introducing an /i-dependent rescaling fac- 
tor, the disk partition function and bulk magnetization can be expressed as 
functions of an effective boundary length and bulk area with no further depen- 
dence on h, except that the bulk magnetization is discontinuous and vanishes 
at h = 0. These results suggest that just as in flat space, the boundary field 
generates a renormalization group flow towards h = oo. An exact analytic ex- 
pression for the boundary magnetization as a function of h is linear near h = 0, 
leading to a finite nonzero magnetic susceptibility at the critical temperature. 



1 



The Ising model with a boundary magnetic field [|I[] has been of renewed interest recently 
as a simple example of a two-dimensional integrable field theory with nontrivial boundary 
interactions . The only boundary conditions for the Ising model which preserve confor- 
mal invariance [f| are free boundary conditions (where the boundary field h vanishes), and 
fixed spin boundary conditions (where h = ±00). Putting an arbitrary field h on the bound- 
ary generates a renormalization group (RG) flow which goes away from the free boundary 
condition towards the fixed boundary condition ||. 

Another subject of recent interest has been the effect of different boundary conditions in 
string theory 0. Just as there are two types of conformally invariant boundary conditions for 
the Ising theory, a conformal field theory of a single bosonic field can have two conformally 
invariant boundary conditions: Neumann and Dirichlet. By considering the continuum limit 
of the Ising model as a single free fermion, in the context of super conformal field theory it 
can be shown that boundary conditions in these two models are related by supersymmetry, 
with Neumann corresponding to free Ising spins and Dirichlet corresponding to fixed spins. 

In this letter we consider the effect of a boundary magnetic field on the Ising model on 
a random surface (the noncritical string with c = 1/2). This theory describes a single Ising 
spin (or equivalently a free fermion) coupled to 2D quantum gravity. The Ising model on 
a random surface can be studied in several ways. One approach is to use the continuum 
formulation of noncritical string theory as a conformal field theory coupled to a Liouville 
field dHH]. Another approach is to describe the model as a matrix model, involving a sum 
over discrete surfaces flO| ; in the discrete formalism, a continuum limit can be taken by 
tuning coupling constants until the surfaces become arbitrarily large; in this limit the theory 
corresponds to the continuous Liouville theory. In this work we will use a discrete formulation 
of the Ising model on a random surface; to calculate correlation functions in this theory we 
use the method of discrete loop equations developed in two previous papers pT|,[12|. Similar 



methods were discussed in [13,14] 



We present here only the results of our investigation. The details of the calculations, 
which are algebraically tedious, will be presented in a later publication. 

In the discrete formulation, the Ising model on a random surface of disk topology has 
a partition function which is given by a sum over all possible triangulations of the disk. 
For each triangulation, the Boltzmann weight is given by placing a single Ising spin on 
each triangle, and summing over all possible spin configurations, giving the Ising partition 



function on that particular geometry. This model can be written as a matrix model |10 

Z(jg,c) = J DUDV exp(-NS(U,V)) , 



with 



S(U, V) = Tr 



-(U 2 + V 2 )-cUV-^-(U 3 + V 3 ) 
2 3 



x 



(2) 
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where U and V are N x N hermitian matrices. 

The first calculation we wish to consider is that of the disk amplitude when the spins on 
the boundary are subjected to an external magnetic field h. In the matrix model language, 
dropping all factors of N as we work in the large N limit, we wish to compute 

oo oo 

4>(h, x, g,c) = J2 ^'(Tr (e h U + e~ h V) k ) = £ x k p k , (3) 

k=0 k=0 

where p k = (Tr (e h U + e V) k ) is the amplitude for a disk with k boundary spins subject to 
the boundary field h. A method for calculating such amplitudes was described in This 
yields a quartic equation satisfied by (ft, in which the coefficients are functions of g, c, h, x, 
and Pi with i < A. We omit this equation for space considerations. 

The quartic equation for <fi gives an exact algebraic solution for the disk partition function 
of the discrete theory. To find the solution in the continuum limit, we must find the critical 
values for x and g at which <ft approaches a singular point. The Ising model is critical for 



the coupling c c = (—1 + 2\/7)/27. The critical value for g is known to be [|K| g c = ylOc^. 
After an analysis of the critical behavior of 0, we find 

x = 9c(l + 2V7)e h 

c e~ 2h + (V7 - 1) + e 2h {2 + V7) ' KJ 

This expression is only valid for h > 0; x c is nonanalytic at h — 0. Throughout the remainder 
of this letter we will restrict attention to the case h > 0; related expressions arise when h < 0. 

To take the continuum limit, we expand around the critical values g = g c e~ e2t , x = x c e~ 6Z . 
Expanding in e, we find 

4/3 

= M*, t, e) + 5 . 27 / Mh f ( Z > T ) + <^ 5/3 ) , (5) 

where <p a is analytic in e. The second term is nonanalytic and describes the behavior in the 
continuum limit. Here 

4/3 / _ r— — x 4/3 



$ = (z+ Vz 2 -AT) + (z - Vz 2 - AT) , (6) 



t is rescaled by a constant factor t = T/5, and z is rescaled by an /i-dependent factor, 
z = a(h)Z where for h > 



e 2M 



a(h) = K —= — '- pr- . (7) 

V ' e~ 2h + (-1 + V7) + e 2h {2 + y/7) y J 

At h = 0, the scaling factor a is discontinuous and goes to a(0) = l/(v / 2 + VlA); the 
constant factor in the universal term in (Kf) also changes discontinuously at this point. Note 
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that the specific form (|7|) for a depends upon the discretization we have chosen for random 
surfaces. 

The universal term in (|5|) can be converted into the asymptotic form of the disk amplitude 
4>(l, a) for fixed boundary length / and disk area a. These forms of the amplitude are related 
through a Laplace transform 

0(4/3) CM) = Jdljda e- zl - ta 4>(l,a) . (8) 
Inverting the Laplace transform, we have 

= _J_LV3 A -7/3 e -^M (9) 

with the rescalings L = a(h)l, A = a/5. Up to an irrelevant multiplicative constant, this 
is precisely the form of the disk amplitude when the boundary conditions are conformal 
|T^ . |T3| . |iI|[T2|j (i.e., with h = or h = ±oo); however, the boundary length I is rescaled by 
the factor a(h) which depends on the boundary magnetic field. 

The boundary magnetization for a spin on the boundary of a disk with k boundary edges 
and n triangles is given by 

, _ (Ti{e h U -e- h V){e h U + e- h V) k - 1 ) n 

{m) ~ (Tr(e h U + e- h V) k ) n ' 1 ' 

where by () n we indicate a sum over triangulations restricted to geometries with n spins (the 
coefficient of g n in an expansion in g). We therefore look at the expectation value of the spin 
at a marked point on the boundary, that is 

oo 

if; = x k+1 (Tr(e h U - e~ h V)(e h U + e~ h V) k ) . (11) 

k=0 

When h = 0, ip vanishes by symmetry. When h ^ 0, we can compute ip(h) by the method of 
loop equations, giving an equation relating ip(h) to 4>{h). Solving this equation, we can find 
the critical expansion of ip{h) about the critical point and the inverse Laplace transform ip 
of the universal part of ip. The boundary magnetization in the continuum limit is then given 
by (for h > 0) 

I = (e h - e- h )(e- h + (2 + V7)e h ) 
{m) (e- 2ft + (-l + v / 7) + (2 + v / 7)e 2ft ) ' 1 J 

Note that the boundary magnetization is independent of / and a. 

A graph of the boundary magnetization is shown in Fig. 1 (bold curve). As expected, 
with no field the magnetization is zero, and for an infinite field the magnetization is 1. This 
result is compared with the boundary magnetization in flat space, computed by McCoy and 
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Wu jlj (dashed curve). Whereas in flat space the magnetization scales as hlnh for small 
h, leading to a divergence in the magnetic susceptibility at the critical temperature, on a 
random surface we find that the magnetization is linear at h — 0, with a finite susceptibility 



X = d h (m)\ h=0 



1 + 2^/7 



(13) 



The two point boundary magnetization can be computed in a similar way, and is found to 
be equal to the square of the one point magnetization. 

<m> 




FIG. 1. Boundary magnetization (m) as a function of boundary field h in flat space (dotted line) 
and on a random surface (bold line) 

Consider now the average bulk magnetization with a boundary magnetic field, on a disk 
with boundary length k and area n. 



<M> 



(Tr(e h U + e~ h V) k Tr(U - V)) n 
(Tr(e h U + e- h V) k Tr(U + V)) n 



(14) 



This can be evaluated by considering cylinder amplitudes with one boundary having a bound- 
ary magnetic field, and the other with a single boundary edge. The second boundary repre- 
sents a marked point on the bulk. Again such a quantity can be computed by the method 



of loop equations ||12|| . After some algebra, it can be shown that the bulk magnetization in 
the continuum limit is given by the simple expression 



(M) = L^A- 1 ^ 



(15) 



Since I and a are measured in lattice units, in the continuum limit L <C A so the magnetiza- 
tion is always less than 1. More precisely, in the continuum limit, (M) scales as <5 1//3 where 
5 is the lattice spacing; this agrees with the scaling dimension of the gravitationally dressed 
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spin operator . Note that this form of the magnetization is independent of h except for 
the dependence through the scaling factor a incorporated in L. At h = 0, this magnetization 
is discontinuous and vanishes. 

We have found that both the disk partition function and the bulk magnetization are 
naturally expressed in terms of a rescaled boundary length L = a(h)l. An interesting 
feature of the bulk magnetization is that, with the particular choice of discretization 
scheme we have used here, when expressed in terms of the actual boundary length I, the 
magnetization is a function which for fixed values of I and a decreases as the boundary 
magnetic field h increases. This counterintuitive result can be explained in terms of the 
influence of the boundary field on the disk geometry. In the vicinity of the disk boundary, 
the existence of a boundary field causes a local fluctuation of the discrete geometry which 
depends upon the magnitude of the boundary field. In the continuum limit, this effect is 
restricted to a vanishingly small region near the boundary. The significance of the rescaled 
boundary length L is that the effects of the boundary magnetic field can be described by 
using in place of the boundary length I an effective boundary length L = a(h)l (see Fig. 2). 
In the bulk of the disk in the continuum limit, for any nonzero boundary field h, all the 
physics is identical to the physics which would occur on a disk of boundary length L/a^ 
with infinite magnetic field on the boundary, where = 1/(2 + v7) is the limit of the 
scaling factor as h — > oo. Since a(h) is a decreasing function of h, the effective boundary 
length L decreases for fixed I as h increases. As in Fig. 2, a decrease in L for fixed A forces 
the disk to deform so that the spins move further away from the boundary, causing a net 
decrease in the average magnetization. 

L 1 

(a) (b) 

FIG. 2. (a) A portion of a surface with area A and boundary length L, (b) A disk with boundary 
length I and effective boundary length L < I is stretched so that an average point is further from 
the boundary 

Another interesting feature of the bulk magnetization is that when expressed in terms 
of the rescaled boundary length L it is independent of h, except where h = 0, when the 
magnetization vanishes. This result indicates that just as in flat space, any nonzero magnetic 
field produces a renormalization group flow whose fixed point limit is the infinite magnetic 
field boundary condition. 
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Finally, it should be noted that of the results presented here, the scaling factor a(h) in 
(0) and the explicit form of the magnetization (IT3) have a functional dependence on h which 
depends on the set of triangulations we have used for the disk. Just as in the case of the flat 
space Ising model, different choices of lattice discretization will give rise to different critical 
values g c , different rescaling functions a(h), and boundary magnetizations with different 
functional dependence on h. Certain of the results obtained here - the functional dependences 
of the disk partition function and the bulk magnetization on the rescaled boundary length 
L and area A, and the fact that the boundary magnetic susceptibility is finite and nonzero 
- should be independent of the choice of discretization, however, and should give the same 
continuum limit in any formulation of the theory. On the other hand, the result that the 
bulk magnetization decreases with increasing boundary field is not necessarily a universal 
result, since it depends upon the explicit form of the rescaling function a(h). 
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